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Abstract 

In our previous article [4] an approach to derive Papapetrou equations for 
constrained electromagnetic field was demonstrated by use of field variational 
principles. The aim of current work is to present more universal technique 
of deduction of the equations which could be applied to another types of 
non-scalar fields. It is based on Noether theorem formulated in terms of 
Cartan' formalism of orthonormal frames. Under infinitesimal coordinate 
transformation the one leads to equation which includes volume force of spin- 
gravitational interaction. Papapetrou equation for vector of propagation of 
the wave is derived on base of the equation. Such manner of deduction 
allows to formulate more accurately the constraints and clarify equations for 
the potential and for spin. 

1 Introduction 

Complete description of motion of a non scalar wave in gravitational field 
is given by its covariant field equations. Under quasi-classical considera- 
tion when wave length is sufficiently less than typical scales of observations, 
propagation of the wave is substituted by motion of a particle. A. Papa- 
petrou in early 50s showed that spin-gravitational interaction changes form 
of the trajectory of a particle with spin. Deduction equations of motion of 
the quantum particle with account of spin-gravitational interaction can be 
provided by construction its classical relativistic Lagrangian. This problem 
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has chained a great interest of researchers for last decades. Nevertheless a 
satisfactory description was not obtained, see [TJ [3] . 

Evident progress in theoretical treatment of the problem for photons was 
achieved in series of our recent articles [2j [3j HJ [5j. Particularly in the 
work [4] an approach to deduction Papapetrou equation for photon based 
on field variational principle was demonstrated. Essence of that approach 
was attachment the field of potential to congruence of 0-curves which spec- 
ify propagation of the wave. The congruence and the field were described 
by frame {n±, ni^} which is orthonormal in sense that < n + , n_ >= 1, 
< n a , n a >= —1, a = 1,2 are only non zero scalar products of the vectors. 
Vector n_ is tangent to the 0-curves everywhere and vectors constitute 
polarization tangent subspace [I]. Integrand of the action was presented as 
quadric form over derivatives of the potential. Under the constraints imposed 
the one was reduced to form similar to action of classic particle with an ad- 
ditional term describing spin-gravitational interaction. Next it was shown 
that varying this action under infinitesimal dragging changing shape of the 
0-curves leads to Papapetrou equation for shape of the 0-curves. Variation 
components of the potential attached to the curves reduced equation for po- 
tential to D_ A a = 0. That in turn led to equation for spin: Z}_ S a p = 0, see 

m 

In this article we are returning to consideration of the problem pursu- 
ing the aim to create more formal and universal approach to deduct quasi- 
classical equations for motion of photon under account of spin-gravitational 
interaction. At the beginning we explain technic of variation in framework of 
theory of exterior differential forms. Next we develop and improve method 
demonstrated in paper [1] . In particular, conception that elements of the or- 
thonormal frame are variables of Lagrangian describing congruence of curves 
of propagation of the wave fruited one of main ideas of the current work. 
That is to consider elements of the frame as variables of gravitational field. 
Then to apply infinitesimal coordinate transformations which transform or- 
thonormal frames, connection coefficients, potential of electromagnetic field 
but leave invariant gravitational field. As usual Noether theorem gives differ- 
ential equations for the field observables like stress-energy tensor (SET) and 
current of spin. However, this time the equations do not exhibit covariant 
conservation law for the SET, but express equation of dynamic of electro- 
magnetic wave. Left hand side (LHS) of the one is usual covariant divergence 
of the SET, but in the right hand side (RHS) instead zero a term interpreted 
as volume force of spin gravitational interaction arises. It happens due to 
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vector character of electromagnetic field. 

Then we develop our quasi-classical approach. It is based on assumption 
of existance a class of locally monochromatic (LM) electromagnetic waves 
which admit geometric optical description. It means that such the wave 
propagates along congruency of worldlines satisfying an equation whose LHS 
coincides with equation of geodesic line in a canonical comoving frame [I]. 
The equation should be derived from the equation of dynamic. Conditions 
which reduce the one to equation for shape of worldlines specify the class of 
LM electromagnetic waves. Analysis LHS of the equation of dynamic gives 
conditions as follows. T ++ should be only non vanishing component of SET of 
the wave and equation of continuity for probability current had to be satisfied. 
In turn RHS of the equation of dynamic should include characteristics of the 
particle but not the ones of the field. It suggests that current of spin can 
be expressed via tensor of spin and vector of propagation: S a b( c ) = $t Sab, 
besides S12 only non zero component of the tensor of spin due to definite 
helicity of photon. By other words we demand that LM electromagnetic 
wave should be transversal wave. Such the requirements are provided by 
transversality of the potential and choice of suitable form of the Lagrangian. 

Analysis of form of the SET and restrictions put onto it gives constraints 
which should be imposed to the potential of the field. The constraints coin- 
cide with the ones were imposed in [4]. Condition of continuity of probability 
current finds its application under deduction of equation for components of 
the potential. Now D_ A a 7^ but proportional to A a , a = 1, 2 and shows 
attenuation of the wave caused by divergence of the beam of electromagnetic 
wave. Nevertheless it does not change our main conclusions obtained earlier 
[HE]. Specifically, under redefinition the spin as normalized V + component 
of current of spin equation for spin leaves the same. 

2 Lagrange formalism in framework of orthonor- 
mal frame 

Action of electromagnetic field in curved space-time usually expressed by 
integral of Lagrangian density £ multiplied to 4-form e of unit volume. For 
aims of our studies it is convenient to express the integrand as whole 4-form 
A which we call 4-form of Lagrangian: 



A 




A 



e = 




(1) 
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where {x 1 } some coordinates in considered domain of space-time and {v a } 
is a field of orthonormal frames dual to {n b }. In general the 4- form of 
Lagrangian: 

A = A(A b ,DA c ,u a ), (2) 

depend on components of field potential A b (x l ), their derivatives and ele- 
ments of orthonormal frame {i> b } which describes gravitational field. Besides, 
it is convenient to represent derivatives of the potential via covariant exterior 
derivatives (CED): 

DA b = (D a A h )v a = dA b - u b c A c , (3) 

where symbol d stands for exterior derivative and uj b c is 1-form of Cartan' 
rotation coefficients. Notion of CED and rules of operating with it should 
be briefly observed. We define CED of any exterior differential form with 
(exterior) tensor indexes as follows: 

i) {r. ,,.:■■''■■■) <ri: jr ■■■>-■ .... 

~< A T„ cJ - b - + ...-uj b c AT fl ...- c - + . . . • 

Usual Leibnitz rules can be easily generalized for the CED under account of 
antisymmetry of exterior product of the differential forms. We also define 
CED of vector: 

Dn a = n b ®u a b - De = n a ®D h {e a )v h , 

where symbol ® stands for tensor product of elements of tangent and cotan- 
gent spaces. Thanks to definitions it is seen that rules as follows are valid: 

n a ®D(T a ) = D{n a ) AT a + n a ®dT a . (4) 

In terms of 4-form of the Lagrangian the field equations have a form as 
follows: 

where d A/d A b is a partial derivative of 4-form of the Lagrangian over field 
variable A b (x). An algebraic derivative d A/d DA b of 4-form A over 1-form 
DA b [6] reduces it to 3-form. CED of the algebraic derivative restores degree 
of the differential form up to 4th. Notion of algebraic derivative provides 
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compact view of formulae especially under consideration variations of the 
4-form of Lagrangian like follows: 

A(. . . , DA + SDA, ...)- A(.... /;.!....) 6 DA A = 5DA A 

5„A = T afe < 5u a , u b >e = T ab 5v a A V = 5z/ a A (6) 

where 5.DA is variation of CED of potential of the electromagnetic field, 
bv a expresses variation of gravitational field variables, T ab is stress-energy 
tensor (SET) of electromagnetic field, asterisk stands for Hodge conjugation 
of exterior differential forms. 



3 Noether theorem and Spin current 

Under infinitesimal rotations of field of the orthonormal frames {v b }'- 

5u a = -e b a u\ e ab + e ba = 0, (7) 
variation of CED of the potential becomes: 

d DA 

SDA c = D5A c + 5uj ab - c -. (8) 

OUJ ab 

Due to field equations, only part of variation of the Lagrangian contributing 
to variation of the action is: 

. dDA c dA 1. dDA c dA 1 . t c 

S^- A — — = -8uj ab - c - A — — = -5uj ab A S ah c V c , 9 

duj ab dDA c 2 duj[ ab] dDA c 2 

where square brackets means antisymmetrization and S ab c *v c stands for cur- 
rent of spin of the field. Variation of the connection 1-form under the con- 
sidered rotations is: 

5u ab = -De ab , 

where De ab = de ab —uj a c e cb —uj b c e ac is CED of tensor of rotations. Substituting 
it into Qj represents variation of as follows: 

5 A = Jd[...] + J l - £ab D[S ab c *v c ] = 0. 

This gives us covariant equation of continuity of 3-form of current of spin: 

D\S ab c *u c ] =0. (10) 
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4 Coordinate variation 



Infinitesimal vector field e drags coordinate hyper- surfaces {x 1 } onto {y 1 } and 
elements of orthonormal covector and vector frames {n a }, {u b } onto dragged 
frames {'n a }, {V 6 }. Applying this variation to action of Gravitational field 
gives equation of covariant continuity for Einstein tensor in natural vector 
frame and in orthonormal one: 

G i ? ]j = D h (G a h .) = 0. 

The same procedure can be performed for ([1]) as follows: 

8A= f A('n 0) %, 'n 3 )'is° A . . .V 3 - / A( n 0) . . >° A...v\ (11) 



where A(no, • • •) means value of 4-form A on 4-vector [ no n\ ri2 ^3] • Under 
this definition variations of variables of the Lagrangian are given by Lie 
derivatives over vector field e. Use of Cartan' first and second structure 
equations expresses variations of v a and u a b as follows: 

5u a = dv a [e ) + de a = De a - u b a { e )u b , (12) 
8oo a b = du a b (e) + d[u b {n ± k 6 A<](e) = Q b {e) + D[u b {e)]\ (13) 

where Q a b = l/2R a b cd u c A v d is 2-form of curvature. This way variation of 
4-form of Lagrangian (j2j) becomes: 

r * p a 9 A A dA . „ 9A 

S A A + 5,A + 5 DA A = 5A b —— + 5DA b A — — + A 7—; 

Noting that variation of CED as before is given by (|5J) where this time 5a; is 
given by (|T3|) we after applying ((Sj) and (J9j) rewrite the above expression as 
follows: 

~ {n ab (e ) + D[u ab ( e )]} A S ab c *v c + D(e a ) A T ab *u b , (14) 



where only first term of variation of v a ([12]) contributes owing to symmetry of 
the SET. Substituting the above expression into ffTTj) we integrate by parts. 
Referring to ffTUj) we obtain: 



Jd[...] + J \l/2Vl ab {e) A S ab c *v c — e a D 



T ab *u b 



0. 
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It brings us equation as follows: 



n c (g> D 



v 



1/2 n c ® R ab cd S ab d e 



(15) 



The one in contrast with analogous equations for scalar fields is not equation 
of covariant conservation of the SET due to non zero RHS appeared. The 
RHS expresses volume force of spin-gravitational interaction while the LHS 
exhibits CED of current of stress-energy-momentum. Because of this we 
will call (fT5l) as equation of dynamic of the electromagnetic wave in curved 
space-time. Accordingly to rule (jlj) LHS of (TT5]) can be represented as: 



fL ® D 



V 



D(n r 



A T c h *v h + n c 



d 



T 



■ i, 



v 



(16) 



It should be marked that use of the above mentioned rule becomes possible 
owing to fact that all indexes in the above expression are dummy. Such form 
of presentation of covariant divergence of the SET is convenient for transition 
to quasi classical consideration. 



5 Quasi classical approach, Schweber Lagrangian 
and constraints 

Conception of our quasi classical description is a deduction of equation for 
vector of propagation n_ from the equation of dynamic (fT5l) . Besides we 
demand that obtained equation becomes coincide with equation of geodesic 
line under switch off the spin. This means that covariant divergence of SET 
should be reduced to covariant derivative of n_. Equality (|T6|) shows that 
this result is provided by constraints as follows: 

T ++ is only non vanishing component of the SET, (17) 



d 



0. (18) 



It should be distinguished that (I17II18P characterize subspace of locally monochro- 
matic (LM) electromagnetic waves and class of canonical comoving frames 
[1] which allows most convenient form of quasi-classical equations for the 
waves. The same form of SET has monochromatic electromagnetic wave in 
flat space-time propagating along n_: 

T ab = 5 a 5 b u 2 \A\\ 
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where u is frequency of the wave. Thus (TT5T) expresses continuity of current 
of probability in case of flat space-time. Under this conditions LHS of (1151) 
becomes: 

D[n-\ AT7 + V" = T~ + D-[nJ\ e. 
This way we reduce ffl5l) to sought form: 

T ++ D_[n_]= n c Rl ab S ahd . (19) 

LHS of the equation coincides with geodesic equation for vector field n_, 
however in the RHS stands a term including contraction of the curvature 
with current of spin of the field which we call force of spin-gravitational 
interaction. Besides the force still is written in terms of field theory. To pass 
to complete quasi classical description we need to determine field Lagrangian. 

Under study process of propagation of LM electromagnetic wave it is rel- 
evant admit Schweber Lagrangian [7] which provides adequate quasi-classical 
description of LM wave: 

C = 1/2 D a A c D b A c < u a , v b > . (20) 

In terms of orthonormal frames we write 4-form of the Lagrangian as follows: 

A = 1/2 DA C ( n a ) {DA C A *u a } . (21) 

Now let's calculate variation of the action under variation of elements of the 
frame: 

5 V K = 1 /2 DA C (5 n a ) DA C A *v a + 1 /2 DA C ( n a ) DA C A 5*u a . (22) 

It is evident that: 5n a = — Su b ( n a ) rib- So first term in ( 1221) can be rewritten 
as follows: 

-1/2 < 5v b ,is d > D b A c D d A c . 
During variation of the second term asterisk conjugation should be varied: 

5 *v a = 1/3! e a bcd 5{u b A v c A v d ) = 1/2 e a bcd 5u b A v c A v d 
v e A 5 *u a = -l/2e a bcd e ecd f 5u b A *v f = {v ae Vbf - S a f 5 e b ) < 5u b , v s > e. 

This brings us: 

<U = - (p a A c D b A c - \j2D e A s D c A^ ab ) < 5u\ v b > e 
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By other words, accordingly to ([61), we find: 

T ab = - (p a A c D b A c - 1/2 D e AfD e A f Vab ) . (23) 

The next part of the variation contributing to RHS of (fT5l) is variation 
over Cartan' rotation 1-form. The one appears in expression of CED of 
potential Q. Evidently b^DAc = —5uj c d A d , so 

5„A = 8uDA c A = 6u c d A A d D b A c *u b = l/25u cd A A [d D b A c] V 6 . 

Now equating this expression to (Q we obtain: 

S cdb = A [d D b A c] . (24) 

Since photons has definite helicity they are presented by electromagnetic 
waves with circular polarization [3]. Hence it is convenient to consider com- 
plex valued amplitudes of the potential which describe phase shift of circu- 
larly polarized waves. For this aim an ordinary procedure of redefinition the 
observables being quadric forms of the potential components is used. That 
is ApA q — > 1/2 A{ p A q j, where curly brackets mean symmetrization over the 
indexes. After applying it to expressions for the SET and the current of spin 
we obtain their form in complex amplitudes: 

T ab = - l - {D a A c D b A c + D a A c D b A c - Vab D c A d D c A d ) , (25) 

= \ (A [d D b A c] + A [d D b A c] ) . (26) 

Now let us find constraints should be put to the potential. As we consider 
LM electromagnetic wave we expect that potential of the wave has a structure 
as follows: 

A = ae i(p , D + = u, £>_0 = O, D + a = (27) 

D + A a = iuA a . (28) 

Due to (!25l) it is seen that condition of dominance of T ++ demands vanishing 
derivatives of the potential along polarization vectors: 

D a A = 0, a = 1,2. (29) 
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Vanishing the derivative along direction of propagation n_ makes the con- 
dition surely satisfied, however it is not only possibility as it will be shown 
later. Anyhow value of T ++ can be calculated now due to fact that r] ++ = 0. 
Under constraints imposed it is reduced to the form of T ++ of monochromatic 
wave in flat space-time: 

T ++ = -uj 2 < A, A >= cu 2 \A\ 2 (30) 

Quasi classical picture assumes that spin of photon moves along field of vec- 
tors n_ and has only S12 component due to fact of definite helicity of photon. 
Constraints (I29II28I) provide validity of the first condition which can be writ- 
ten as S c db = Scdrjb-, where S c d is tensor of spin. The second condition is 
supplied by additional constraints: 

A± = 0. (31) 



6 Field equations 



Under varying Lagrangian (l2"Uj) over components of potential we obtain re- 
duced form of Maxwell equations: 



D 



D h A c V 6 



0. 



Let's consider equations for nonzero components of the potential. After im- 
posing constraints we obtain: 



D {+ (D- } A a )e+J2 D c A a d*u c = 0. 



(32) 



c=± 



As usual under quasi classical (geometric optical) consideration we assume 
that values of second derivatives of amplitudes of the potential should be 
neglected: 



uo- l \D a (D b a c )\ =u- 1 0[\a 



D_ \a\ 



(\a\~ l D^ 



It reduces (|32l) to: 

[2iuD_ + i(D_ u)] A a t + iuoA a d V + + (£>_ A a ) d *v~ = 0. (33) 
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Complex conjugation of (13"3l gives equation for complex amplitude. Let's 
consider contraction of (13"3"j) with A a and contraction of the one complex 
conjugated with A a . Difference between them gives: 



ujD_ < A, A > +(£>- lj) < A, A > 



e + ito < A,A> dV 



+ 1/2 



A n D_ A a - A n D_ A c 



d*W 



0. 



Due to structure of potential and constraints we daresay that last term in 
the above equation is zero. This manner we rewrite the equation as follows: 

uj- l d[to 2 < A, A > V+] - (D- u)<A,A >= 0. 

But owing to fl30l) and (fTS|) first term in the equation is zero. It gives £)_ uj = 
0. By other words value of derivative of the frequency in chosen frame along 
vector n_ is vanishing in quasi classical approximation: 



D_uj = O 



D- < a, a > 
< a, a > 



-i 2^ 



u D_ < a, a > 
< a, a > 



(34) 



This result allows to rewrite fl33l) and ffl8|) as follows: 



2D_ A a e + A a d*v + 



'D_ A a )*d*v~ = 0, 



d 



< A, A > *u + 



0. 



(35) 
(36) 



Equating to zero factor at u in (135!) as it is used to do under geometric 
optical approximation we obtain equation for potential components. In turn 
(|36|) expresses conservation of probability current of the wave and allows to 
exclude d*v + from the equation for the potential: 



D_A a = 1/2 A a *d*u + = A 



2A r A c 



(37) 



Equation (137j) describes attenuation of the amplitude of potential caused by 
divergence of worldlines of propagation of the wave. 
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7 Explicit form of SET, current of spin and 
Papapetrou equation 



Substituting (j37j) together with constraints into f)25|26p allows us to obtain 
an explicit form of elements of SET: 



D a A b D a A l 



D{ + A h D-\A 



-iuj A b A b /2 + l/2A b iuA b \ 
0, T__ 



D_ A r A c 



A C A C 

-D- A b D- A b < T + 



0. 



= -1/2 D {+ A C D_ } A C -u, , - ^, +H 

The calculations confirm assumptions about view of elements of SET of LM 
wave. Next step is to calculate elements of spin tensor. Thanks to constraints 
all its elements with third index different from "±" is surely zero. But 
accordingly to fl26l) : 

S c d+ = i&A[d A c ] —: Sd, 



S c d- 
Scdb 



D-\A\- 
\A\ 2 

Sfj S cd . 



(A [d A c] + A [d A c] ) =0, 



(38) 



It approves expected factorized form of current of spin. Now substituting 
(1301) and (1581 into ( 1T91 we obtain equation for vector n_ being tangent to 
worldlines of propagation of the wave: 

u 2 \A\ 2 D_ n_ = 1/2 n c R abc+ S ab . (39) 

In fact the one coincides with Papapetrou equation for trajectory of photon 
obtained us in work {!], although combinatorial factor 1/2 was not present 
there. It happened only by a variety in the definitions of tensor S ab in the 
articles. Under practical calculations it is more convenient to introduce nor- 
malized tensor of spin [5]: 

Sab _ A[ a A b ] 



Cab 



oj\A\ 



i\A 



(40) 



As expected, the one is transporting parallel itself in polarization tangent 
subspace: 



ai3 



D_ {A [a )A p] + A [a D_ {A ji]) A[ a Ap] D„ (A C A C 



A C A C 



D_ (A C A C ) r , . - A . . 



{A d A d f 
A [a A p] D_ {A C A C 



2(A d A° 



(A d A 



d\2 



0. 
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Finally excluding potentials of the electromagnetic field, we write system of 
Papapetrou equations for vector of propagation of LM electromagnetic wave 
and its spin as follows: 

D_a Q/3 = 0, ujD. (n_) = 1/2 n c R af3 c _a a ^; (41) 

where value of frequency to leaves constant in canonical comoving frame along 
the trajectory with accuracy up to second order in accordance with j3]. 

8 Remarks 

It should be underlined that idea of this work is grounded on conception of 
article [I]. Specifically, in the above cited paper we started by consideration 
gauge invariant Lagrangian C = 1/2 D a A b D c A d < v a f\v b , v c f\v d >. Next we 
expand scalar products as < u a , z/' c >< u b , > and put the constraints. It 
was guessed that only terms < u a , u c >< u b , v d > would contribute to variation 
of the action due to fact that they yields physically interpretable expressions. 
In present work that assumption finds its expression in choice of form of 
Lagrangian. From another point of view selection of Schweber Lagrangian 
together with the constraints serves transversality of LM electromagnetic 
wave. 

We also note that variation of omitted scalar product in Schweber La- 
grangian D a A c DbA c — > D a A c DbAd < z/ c , v d > may contribute only to the 
SET. Calculations with account of the constraints and symmetrization over 
complex amplitudes analogous to the ones performed in section [7] shows that 
the contribution vanishes. 

Equation for probability current (136|) can be rewritten as: 

D- < A, A > + < A, A > Div ri_ = 0, 

where Div stands for covariant divergence of the vector. Under our ap- 
proximation we daresay that both terms in the equation are small but not 
vanishing. But in our former work [4] it was implicitly assumed that (13T)|) is 
satisfied by vanishing of both of its terms. Hence equation for amplitudes of 
the potential obtained us earlier is particular case of (|37|) . The last is realized 
when worldlines of propagation of the wave are parallel themselves locally, 
as it was demanded in the previous article. 
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